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High-quality predictions for neutrino-nucleus 
cross sections are needed for use in on-going ex- 
perimental studies of neutrino oscillations at GeV 
energies, where a fully relativistic treatment of 
the neutrino-nucleus scattering is mandatory, but 
hard to achieve. Of course any reliable calcula- 
tion for neutrino scattering should first be tested 
against electron scattering in the same kinemati- 
cal conditions. In fact, while relativistic modeling 
of the nuclear dynamics is capable of getting the 
basic size and shape of the inclusive cross section, 
it can hardly account for important details of the 
response. Specifically, the data display a long tail 
which is largely absent in most modeling, possibly 
due to the absence of classes of short-range cor- 
relation effects and to the treatment of final-state 
interactions. 

In order to overcome these difficulties, instead 
of using a specific nuclear model, we propose here 
a method for extracting informations on the nu- 
clear dynamics from electron scattering exper- 
iments and use them to predict the neutrino- 
nucleus cross section. 

The approach is based on the superscaling be- 
havior of the electron-nucleus cross section in 
both the quasi-elastic (QE) and A peaks. De- 
tailed studies of superscaling can be found in 
Refs. |1I2| ; here we simply recall the basic idea. 
By dividing the experimental inclusive cross sec- 
tion for various momentum transfers q and nuclei 
(or, equivalently, Fermi momentum fcj?) by an ap- 
propriate single-nucleon cross section, a reduced 
cross section /, called the superscaling function, 
is obtained, which embodies the nuclear dynam- 



ics. This can be plotted versus a well-chosen scal- 
ing variable: if no dependence of / upon q is ob- 
served, we call this behavior scaling of the first 
kind, whereas if no dependence occurs on the spe- 
cific nucleus, we call it scaling of the second kind; 
if both types of scaling behavior are found we say 
that superscaling occurs. 

Several choices for the appropriate scaling vari- 
able have been proposed in the literature. Here 
we choose that which naturally emerges from 
the relativistic Fermi gas (RFG) model. In this 
framework the nuclear response functions in both 
the quasi-elastic (X = QE) and A-resonance 
(X = A) regions have the general structure 



Rx 



NrriN 
qk F 



RxY 



(1) 



where M is the appropriate nucleon number and 
[RxY is the single-nucleon response function. 
The RFG superscaling function 



■tf)9(l-ip) (2) 



is independent of q and hp, being a function of a 
specific combination of q, lu and kp , the scaling 
variable, defined as 
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where Tp is the Fermi kinetic energy, Arnj^T = 
zndpx = l + (m x -m%)/(q 2 -^) (m QE 
being the nucleon mass). The sign in Eq. © is 
+ (— ) if the energy transfer to is higher (lower) 
than the one corresponding to the maximum of 
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the peak. Moreover, a small energy shift is in- 
troduced, implying ip — > ip', in order to force the 
maximum of the response to occur for ip' x = 0. 

At sufficiently high energies the electron scat- 
tering data exhibit both types of scaling behavior 
in the quasi-elastic region. For specific nuclei the 
first-kind scaling is quite good at excitation en- 
ergies below the QE peak (the so-called scaling 
region): this is the familiar y-scaling behavior. 
At energies above the peak, where nucleon reso- 
nances (especially the A) are important, this type 
of scaling is broken. The scaling violations appar- 
ently reside in the transverse response, but not in 
the longitudinal, which appears to superscale. In 
fact, this is not unexpected, since MEC and in- 
elastic contributions, which do not scale, are pre- 
dominantly transverse in the kinematic regions of 
interest. On the other hand scaling of the second 
kind works very well in the scaling region and, 
even in the resonance region, is only violated at 
roughly the 20% level. 

The approach taken in |5] therefore has been 
to use the experimental longitudinal responses to 
define the superscaling function. Reliable sepa- 
rations of data into their longitudinal and trans- 
verse contributions for A > 4 are available only 
for a few nuclei; all of these response functions 
have been used to extract the "universal" quasi- 
elastic response function fQ E and to obtain a pa- 
rameterization by a simple function. Fig. 1 (up- 
per panel) shows f® E (ipQ E ) averaged over the 
nuclei employed, together with the corresponding 
fit, plotted against the scaling variable iP'qe- 

Note that / < ^ e (V'q_b) nas a somewhat asym- 
metric shape and a tail that extends towards 
positive values of iP'qe- In contrast, the RFG 
superscaling function J2J is symmetric, is lim- 
ited strictly to the region — 1 < i/j'qe < +1 and 
has a maximum value of 3/4, while the empiri- 
cal superscaling function reaches only about 0.6. 
The microscopical origin of this asymmetric shape 
has been the object of recent investigations in 
the context of the relativistic impulse approxi- 
mation Let us now focus on the resonance 
region. In order to isolate the A contribution, we 
subtract from the total experimental cross sec- 
tions the quasi-elastic cross section recalculated 
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Figure 1. Experimental QE and A superscaling 
functions versus the corresponding scaling vari- 
ables 

using the universal f < * E {ip , QE) introduced above. 
That is, we remove the impulsive longitudinal and 
transverse contributions that arise from elastic 
eN scattering, leaving (at least) MEC effects with 
their associated correlations and impulsive contri- 
butions arising from inelastic eN scattering. The 
MEC effects can be ignored as they provide rela- 
tively small corrections 0j , and thus this yields a 
response which is largely dominated by the A, at 
least for energy losses lower than the maximum 
of the A peak. 

Using the same procedure as in the QE peak, 
we can again reduce the left-over cross section by 
dividing by the appropriate single-nucleon N — > 
A cross section and display the result versus the 
scaling variable ip' A defined in Eq. ©. In Fig. 1 
(lower panel) we show the resulting / a (V , a) ex ~ 
tracted from the high-quality world data for in- 
clusive electron scattering from 12 C and 16 O in 
the QE and A regions. These data span energies 
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extending from 300 MeV to 4 GeV and scattering 
angles from 12 to 145 degrees, depending on the 
beam energy. As for f^ E , the experimental val- 
ues of f A have been parameterized by a simple 
analytical function (solid line). The data appear 
to scale reasonably well up to the peak of the A 
(ip' A = 0), although clearly for still higher exci- 
tation energies the scaling is broken by processes 
that are not well represented via A-dominance 
(other resonances and, at the larger values of q, 
the tail of deep inelastic scattering). There is also 
some excess at large negative ip' A which breaks 
the scaling to some degree, probably due to con- 
tributions from MEC and their associated corre- 
lations @|. 

The above analysis can now be checked by as- 
sembling all of the pieces obtained via the scal- 
ing procedures to produce a total inclusive cross 
section which can be compared with data. Sev- 
eral examples of such a comparison are shown in 
Rcf. 5]. The answers are very encouraging, and 
only for specific kinematics do we see deviations 
as large as 10-20%, probably due to MEC and 
correlations. In particular the dip region, usually 
hard to get from microscopic models, is nicely fit- 
ted in the present context. 

We are then in a position to take a step in a 
new direction. Since the above superscaling func- 
tions, upon being multiplied by the electromag- 
netic N — > N and N — > A single-nucleon cross 
sections, reproduce quite well the total nuclear 
electromagnetic cross section, we can just as well 
multiply by the corresponding charge-changing 
weak interaction N — > N and N — > A single- 
nucleon cross sections to obtain predictions for 
neutrino reactions in nuclei at similar kinemat- 
ics. The result for the reaction 12 C(i/ M ,/i~) is 
shown in Fig. 2 versus the muon momentum k' for 
typical kinematical conditions. Cross sections at 
different kinematics and for antineutrino scatter- 
ing can be found in Ref. 0. As discussed above, 
the predictions at momenta to the left of the A 
peak (excitations lying above the A region) are 
unreliable, since our scaling approach does not 
fully account for meson production, including res- 
onances other than the A, and DIS processes. For 
comparison, in Fig. 2 we also show the cross sec- 
tion obtained using the RFC Modeling based on 
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Figure 2. v- 12 C cross section plotted versus the 
muon momentum k! . Dash-dotted: QE; dashed: 
A; solid: total; heavier lines (red on-line): RFC 

mean-field theory where finals state interactions 
are not included is very similar to the RFG result, 
since at the relatively high energies of interest the 
dynamical effects embodied in an effective mass 
are expected to be small. Likewise, relativized 
shell model predictions are close to the RFG pre- 
dictions [d] and, moreover, RPA correlations are 
expected to be relatively small for the high ener- 
gies involved. Thus, the RFG predictions effec- 
tively represent a larger set of models. As can be 
seen in the figure, all therefore differ significantly 
from the scaling predictions. 

In conclusion, given the success of the scal- 
ing approach in studies of inclusive electron scat- 
tering for the kinematic region under study, we 
expect neutrino reaction cross sections also ob- 
tained using scaling ideas to be more robust than 
those based directly on existing models. 
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